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Abstract 
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O Hadron generation in the Geant4 simulation tool kit is compared with inclusive spectra 

^ of secondary protons and pions from the interactions with beryllium nuclei of +8.9 GeV/c 

^ protons and pions, and of —8.0 GeV/c pions. The data were taken in 2002 at the CERN 

^ Proton Synchrotron with the HARP spectrometer We report on significant disagreements 

between data and simulated data especially in the polar-angle distributions of secondary 

protons and pions. 
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1 Introduction and motivation 

The HARP experiment arose from the realization that the differential cross- sections of hadron 
production in the collisions of few GeV/c protons with nuclei were known only within a factor of 
two to three. Consequently, the HARP spectrometer was designed to carry out a programme of 
systematic and precise measurements of hadron production by protons and pions with momenta 
from 3 to 15 GeV/c. The experiment was in operation at the CERN Proton Synchrotron in 2001 
and 2002, with a set of stationary targets ranging from hydrogen to lead, including beryllium. 

The data from the HARP spectrometer can be used, amongst other purposes, for the physics 
validation of hadron generators that are used in simulation tool kits such as Geant4 [1]. This is 
of interest for the correct interpretation of data that will be forthcoming, e.g., from experiments 
attheLHC[2]. 

In this paper, data are used from the HARP large-angle spectrometer that comprised a cylin- 
drical Time Projection Chamber (TPC) and an array of Resistive Plate Chambers (RPCs) around 
the TPC. The purpose of the TPC was the measurement of the transverse momentum p^ and of 
the polar angle 9 of tracks, and particle identification by dE/dx. The purpose of the RPCs was 
a complementary particle identification by time of flight. 

The data analysis that underlies the spectra shown in this paper rests on the calibrations of 
the TPC and the RPCs that our group published in Refs. [3] and [4]. For a more detailed account 
of our calibration work we refer to our collection of memos and analysis notes [5]. We recall 
that we disagree with the calibrations and physics results reported by the 'HARP Collaboration' , 
as discussed in Refs. [6] and [7]. 

With a view to correcting for losses of secondary particles from acceptance cuts, and for 
migration due to finite detector resolution, the measurement of tracks in the detector must be 
simulated with a Monte Carlo program. We use the Geant4 tool kit for this purpose. It was at 
this point that we noticed peculiar structures in the polar-angle spectra of secondary particles 
generated by Geant4's LHEP 'physics list' that prevented the weighting of generated tracks 
by smooth functions. Further investigations showed that this is a rather common phenomenon 
across Geant4's hadronic physics lists. 

Figure [T] shows a typical example of an unphysical structure in generated longitudinal mo- 
mentum pl versus transverse momentum p^ of secondaries. Since the structure is genuinely 
connected with the polar angle 9, it tends to be washed out when integrating over either pl 
or Pt^^- That may explain why — as nearly as we can tell — these structures were not noticed 
before. 

2 Hadron generators in the Geant4 simulation tool kit 

The Geant4 simulation tool kit provides several physics models of hadronic interactions of 
hadrons with nuclei, and collections of such models, termed physics lists. The latter are tailored 
with a view to optimizing performance for specific applications. 

Table[T]lists and characterizes a representative selection of physics lists of hadronic interac- 
tions in Geant4^\ together with the used physics models and the energy ranges where the latter 
are considered to be reliable [8]. 

In the so-called 'low-energy' domain (defined as kinetic energy E of the incoming hadron 
below some 25 GeV), a modified version of the GHEISHA package of Geant3 is used in many 
physics hsts: the Parametrized Low-Energy Model ('LE GHEISHA'). Optionally, for E below 



^^ All physical quantities in this paper refer to the laboratory system. 
2) Version 9.1 dated 14 December, 2007. 
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Fig. 1: Longitudinal momentum pL versus transverse momentum p^, as generated by Geant4's LHEP 
physics list for secondary 7r+ from the interactions of +8.9 GeV/c beam tt"*" with beryllium nuclei at 
rest. 



a few GeV, the Bertini Cascade [9] ('BERT') or the Binary Cascade [10] ('BIC') models can be 
enabled, with a view to simulating the cascading of final-state hadrons when they move through 
nuclear matter. As an alternative to LE_GHEISHA, a modified version of the FRITIOF string 
fragmentation model [11] ('FTF') is available. 

In the so-called 'high-energy' domain, mostly the Quark-Gluon String Model ('QGSM') 
is used, with FTF and the Parametrized High-Energy Model ('HE_GHEISHA') as alternatives. 
Further terms that appear in Table [T] and are explained in Ref. [8], are 'PRECO' for the Pre- 
compound model, 'QEL' for the Quasi-elastic scattering model, and 'CHIPS' for the Chiral 
Invariant Phase Space model. 

The energy ranges of models tend to overlap. In the overlap region, the model is chosen 
randomly but the choice is biased by the difference between the kinetic energy of the beam 
particle and the kinetic energy limits of the models. 

Below, we compare the predictions of Geant4 hadronic physics lists with our data: the 
inclusive proton, 7r+ and vr^ spectra that are generated by the interactions with beryllium nuclei 
of +8.9 GeV/c protons and 7r+, and of -8.0 GeV/c tt". 



3 The HARP large-angle spectrometer 
3.1 Physics performance 

For the purpose of this paper, the essential physics performance parameters are the resolution 
and the scale of the transverse momentum p^ of final- state particles, the resolution and the scale 
of the polar angle 9, and the separation of pions from protons. We briefly give evidence of the 
salient features, and refer the reader to our respective technical publications [3,4] for details. 

The resolution of the inverse transverse momentum measured by the TPC depends slightly 
on the relative velocity /3 and on 6 of the particles. It is in the range 0.20 < a{l/p'j') < 
0.25 (GeV/cy^. Figure [2] shows the difference of the inverse transverse momentum of positive 



Table 1 : Overview of selected physics lists of hadronic interactions in Geant4. 



Physics Hst 


Proton beam 


vr^ beam 


LHEP 


HE.GHEISHA 
LE_GHEISHA 


25 GeV-100 TeV 
0-55 GeV 


HE.GHEISHA 
LE_GHEISHA 


25 GeV-100 TeV 
0-55 GeV 


LHEP_PRECO_HP 


HE_GHEISHA 

LE.GHEISHA 

PRECO 


25 GeV-100 TeV 
0.15-55 GeV 
0-0.17 GeV 


HE_GHEISHA 
LE.GHEISHA 


25 GeV-100 TeV 
0-55 GeV 


QGSC 


QGSM+QEL+CHIPS 
LE.GHEISHA 


8 GeV-100 TeV 
0-25 GeV 


QGSM+QEL+CHIPS 
LE.GHEISHA 


8 GeV-100 TeV 
0-25 GeV 


QGS_BIC 


QGSM+BIC 

LEGHEISHA 

BIG 


12 GeV-100 TeV 
9.5-25 GeV 
0-9.9 GeV 


QGSM+QEL 

LEGHEISHA 

BIG 


12 GeV-100 TeV 
1.2-25 GeV 
0-1.3 GeV 


QGSP 


QGSM+QEL+PRECO 
LEGHEISHA 


8 GeV-100 TeV 
0-25 GeV 


QGSM+QEL+PRECO 
LEGHEISHA 


8 GeV-100 TeV 
0-25 GeV 


QGSPBERT 


QGSM+QEL+PRECO 

LE.GHEISHA 

BERT 


8 GeV-100 TeV 
9.5-25 GeV 
0-9.9 GeV 


QGSM+QEL+PRECO 

LE.GHEISHA 

BERT 


8 GeV-100 TeV 
9.5-25 GeV 
0-9.9 GeV 


QGSPBIC 


QGSM+QEL+PRECO 

LE.GHEISHA 

BIG 


8 GeV-100 TeV 
9.5-25 GeV 
0-9.9 GeV 


QGSM+QEL+PRECO 
LE.GHEISHA 


8 GeV-100 TeV 
0-25 GeV 


QBBC 


QGSM+QEL+CHIPS 
BIG 


6 GeV-100 TeV 
0-9 GeV 


QGSM+QEL+CHIPS 
BERT 


6 GeV-100 TeV 
0-9 GeV 


FTFC 


FTF+QEL+CHIPS 
LEGHEISHA 


4 GeV-100 TeV 
0-5 GeV 


FTF+QEL+CHIPS 
LEGHEISHA 


4 GeV-100 TeV 
0-5 GeV 


FTFP 


FTF 
LE.GHEISHA 


4 GeV-100 TeV 
0-5 GeV 


FTF+QEL+PRECO 
LE.GHEISHA 


4 GeV-100 TeV 
0-5 GeV 


FTFP_BERT 


FTF 
BERT 


4 GeV-100 TeV 
0-5 GeV 


FTF+QEL+PRECO 
BERT 


4 GeV-100 TeV 
0-5 GeV 



particles with 0.6 < /3 < 0.75 and 45° < 9 < 65° from the measurement in the TPC and 
from the determination from RPC time of flight with the proton-mass hypothesis. The positive 
particles are protons, the background from pions and kaons is very small. Subtracting quadrati- 
cally from the convoluted resolution of 0.27 (GeV/c)^^ the contribution from the time-of-flight 
resolution of the RPC, gives a net TPC resolution of (t(1/]9t) = 0.20 (GeV/c)^^. 
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Fig. 2: Difference of the inverse transverse momenta of positive (shaded histogram) and negative (black 
points) particles from the measurement in the TPC and from the determination from RPC time of flight, 
for 0.6 < j3 < 0.75 and for 45° < 6 < 65°; the positive particles are protons, with a very small 
background from pions and kaons. 



From the requirement that 7r+ and vr with the same RPC time of flight have the same mo- 
mentum, the momentum scale is determined to be correct to better than 2%, for both positively 
and negatively charged particles. 

The polar angle 9 is measured in the TPC with a resolution of ~9 mrad, for a representative 
angle of 6' = 60°. To this a multiple scattering error has to be added which is ~7 mrad for a 
proton with p^ = 500 MeV/c and 9 = 60°, and ~4 mrad for a pion with the same characteristics. 
The polar-angle scale is correct to better than 2 mrad. 

As for the separation of pions from protons: the resolution of dE/dx in the TPC is 16% 
for a track length of 300 mm, and the system time-of-flight resolution is 175 ps. Figure |3] 
(a) shows the specific ionization dE/dx, measured by the TPC, and Fig. [3] (b) the relative 
velocity /? from the RPC time of flight, of positive and negative secondaries, as a function of the 
momentum measured in the TPC. The figures demonstrate that, in general, protons and pions 
are well separated. They also underline the importance of the complementary separation by 
RPC time of flight at large particle momentum. The average values of dE/dx and (3 agree well 
with theoretical expectations, thus confirming the validity of the absolute scales of momentum, 
dE/dx, and time of flight. 
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Fig. 3: Specific ionization AElAx [in units of minimum-ionizing pulse height] (a) and velocity (5 (b) 
versus momentum [GeV/c], for positive and negative tracks in +8.9 GeV/c data. 



3.2 Acceptance and migration 

We discuss differences of inclusive spectra between data and simulated data in terms of the 
distribution in the polar angle 6, for different ranges of px- The primary reason for this choice 
is that disagreements show up most clearly in 6. At the same time, 6' is a well-measured exper- 
imental quantity. We also consider that the 9 distribution provides the clue to the origin of the 
disagreements. 

The physics performance parameters in the transverse momentum px and polar angle 6 
are so good that finite resolution, or a small dependence of acceptance cuts on px or 0, does 
not appreciably affect the comparison of data with simulated data (the chosen ranges of px 
exceed by a factor of two or more the px resolution, and the chosen bin size of 2° (35 mrad) 
of 6 exceeds by a factor of two the 6 resolution). That is substantiated in Fig. |4] which shows 
examples of the dominant disagreements between data and simulated data: patterns reminiscent 
of diffractive scattering [Fig. |4](a)] and of elastic scattering [Fig. |4](b)]. The full lines show 
the Monte Carlo generated polar-angle distributions, while the crosses show the same after 
acceptance cuts, particle identification cuts, and with resolution effects included. Obviously, 
the structures seen in the simulated data are not appreciably altered by experimental acceptance 
and resolution. 

We conclude that for the comparison of the shapes of inclusive particle spectra between data 
and simulated data — which is the purpose of this paper — it is sufficient to compare data with 
Monte Carlo generated distributions, without correction of losses from acceptance cuts and of 
migration stemming from finite resolution. Also, the comparison is intentionally restricted to 
kinematical regions where there is ample and unambiguous separation of pions from protons'^^ 
To identify a secondary particle as a pion or as a proton it is required that the measured AEIdx 
and time of flight are both consistent with the given particle hypothesis. In addition, the time 
of flight has to be inconsistent with the opposite hypothesis. For particles without either AElAx 
or time-of-flight measurement, the cut on the available variable is tightened. Within the ac- 



■^' Absolutely normalized double-differential cross-sections, obtained after due corrections for acceptance and 
migration, and making use of proper weights for particle identification and therefore spanning larger ranges of 
kinematical parameters than discussed in this paper, will be presented in forthcoming papers. 
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Fig. 4: Comparison of Monte Carlo generated (full lines) with Monte Carlo accepted (crosses) tracks; 
(a) polar-angle distributions of vr^ for incoming 7r+, and (b) of protons for incoming tt"*". 



cepted phase space, the particle identification efficiency is between 70% and 90%, while the 
contribution from wrong particle identification is below 5%. 

4 Data versus simulation from selected Geant4 physics lists 

The combination of the choice of hadron generators with the choice of incoming beam particles 
and the choice of secondary hadrons, leads to a large a number of possible plots. With a view 
to simplifying matters, we select for several physics lists two plots each that are representative 
for the agreement and disagreement, respectively, between data and simulation. 

In all plots, we compare the Monte Carlo-generated 9 distribution with the 9 distribution of 
data. Positive beam particles have +8.9 GeV/c momentum, and negative beam particles have 
—8.0 GeV/c momentum. The target is a 5% Aabs thick stationary beryllium target. In Figs. |5] 
to [lOl data are shown as crosses while simulated data are shown as full lines. As justified in 



Section 3.2 the data are not corrected for losses from acceptance cuts and migration stemming 



from finite resolution. With a view to emphasizing shape differences, the data are normalized 
to the simulated data in the angular range 20° < 9 < 125°. 

Figures [5] to [To] show the comparisons for several physics lists from Table [Tj LHEP, QGSC, 
QGSP_BERT, QGSP_BIC, QBBC, and FTFP 

There are three distinct problems visible in the comparison of 9 distributions of data and 
simulated data: 

- an unphysical peak for secondary protons near 9 = 70°; 

- an unphysical diffraction-like pattern for secondary pions; 

- a poor agreement in the shape. 

The different physics lists behave differently with respect to the type of disagreement, yet 
the unphysical peak for secondary protons near 9 = 70°, and the unphysical diffraction-like 
pattern for secondary pions appear as dominant problems. 
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Fig. 5: LHEP physics list; polar-angle distributions of protons for incoming protons (left panel), and of 
7r+ for incoming 7r+ (right panel). 
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Fig. 6: QGSC physics list; polar-angle distributions of protons for incoming protons (left panel), and of 
vr"'" for incoming protons (right panel). 
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Fig. 7: QGSP_BERT physics list; polar-angle distributions of n for incoming 7r+ (left panel), and of 
protons for incoming tt^ (right panel). 
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Fig. 8: QGSP_BIC physics list; polar-angle distributions of 7r+ for incoming protons (left panel), and of 
7r~ for incoming tt"*" (right panel). 
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Fig. 9: QBBC physics list; polar-angle distributions of protons for incoming 7r+ (left panel), and of tt 
for incoming 7r+ (right panel). 
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Fig. 10: FTFP physics list; polar-angle distributions of protons for incoming 7r+ (left panel), and of vr^ 
for incoming tt^ (right panel). 
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With a view to elucidating the physics origin of these dominating problems, we examine 
more closely for the LHEP physics list the disagreements between data and simulated data for 
different combinations of incoming and secondary particle. 



Figure 1 1 shows for a specific range of p^ all combinations of incoming and secondary 
protons, 7r+ and n^ . 
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Fig. 11: LHEP physics list; polar-angle distributions of protons (left panels), 7r+ (middle panels), and 
vr^ (right panels), for incoming protons (top row), incoming 7r+ (middle row), and incoming vr^ (bottom 
row). 



Figure \T2\ presents for the LHEP physics list in four ranges of px the spectra of secondary 



pions from incoming protons, and Fig. 13 the same for secondary protons from incoming pions. 
The conclusions from the comparison of data with simulated data are summarized in Table [2j 
We qualify the various physics lists in the order 'good', 'acceptable', 'poor', 'unacceptable'. 
Note that our qualification refers to the restricted polar-angle range 20° < 9 < 125°. 
None of the standard physics lists of Geant4 is qualified 'good' . 
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Fig. 12: LHEP physics list; polar-angle distributions of Tr"*", for incoming protons, in four different ranges 
of pt- 



The unphysical peak in the distribution of secondary protons around 70° is consistent with 
the kinematics of elastic scattering of the incoming particle with a proton at rest. The diffraction- 
like pattern in the distribution of seconday pions is consistent with diffractive scattering of the 
incoming particle on a stationary disc with the diameter of a nucleon. We conjecture that the dif- 
ferences between data and simulated Geant4 data arise dominantly from an inadequate descrip- 
tion of the elastic scattering and diffraction scattering of incoming beam particles on nucleons 
embedded in a nucleus. 

For the analysis of our data, we have used for incoming beam protons the QGSP_BIC 
physics list and see no strong reason to reconsider this choice. For incoming beam pions, 
none of the standard physics lists for hadronic interactions was acceptable, so we had to build 
our private HARP_CDP physics list. This physics list starts from the QBBC physics list (see 
Table [T]). Yet the Quark-Gluon String Model is replaced by the FRITIOF string fragmentation 
model for kinetic energy E > 6 GeV; for i? < 6 GeV, the Bertini Cascade is used for pions. 
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350 < p-^ < 450 IVIeV/c 




600 



20 40 60 80 100 120 140 160 180 

Polar angle (degrees) 



Pliysics list: LHEP 
:"*" beam 
secondaries p 
450 < p-^ < 550 MeV/c 




20 40 60 80 100 120 140 160 180 

Polar angle (degrees) 



Physics list: LHEP 
Tt"*" beam 
secondaries p 
550 <Pj< 700 MeV/c 




Physics list: LHEP 
n"*" beam 
secondaries p 
700 < p-;- < 850 MeV/c 




20 40 60 80 100 120 140 160 180 

Polar angle (degrees) 



20 40 60 80 100 120 140 160 180 

Polar angle (degrees) 



Fig. 13: LHEP physics list; polar-angle distributions of protons, for incoming vr^, in four different ranges 

of pt- 

and the Binary Cascade for protons; elastic and quasi-elastic scattering is disabled^^ 

We qualify our HARP_CDP physics list as 'acceptable' in all four categories listed in Ta- 
ble [H 

5 Summary 

We have presented significant disagreements in the laboratory polar-angle distributions between 
data from the HARP large-angle spectrometer, and data simulated by various physics lists in 
the Geant4 simulation tool kit. An unphysical peak for secondary protons near 9 = 70°, and an 
unphysical diffraction-like pattern for secondary pions appear as dominant problems. 



^'Our Monte Carlo simulation will be discussed in the necessary detail in a forthcoming paper. 
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Table 2: Conclusions on Geant4 standard physics lists of hadronic interactions. 



Physics list 


Proton beam 


n"^ beam 


Secondary 
protons 


Secondary 

7r± 


Secondary 
protons 


Secondary 


LHEP 


poor 

(shape) 


unacceptable 

(diffr. patt.) 


poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


LHEFPRECOJiP 






poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


QGSC 


poor 

(shape) 


unacceptable 

(diffr. patt.) 


poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


QGS_BIC 






poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


QGSP 


poor 

(shape) 


unacceptable 

(diffr. patt.) 


poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


QGSPJBERT 


poor 

(shape) 


poor 

(shape) 


poor 

(el. scatt. peak 
in 7r+ beam) 


poor 

(shape) 


QGSPJBIC 


poor 

(shape) 


poor 

(shape) 


poor 

(el. scatt. peak 
in 7r+ beam) 


unacceptable 

(diffr. patt.) 


QBBC 


unacceptable 

(el. scatt. peak) 


acceptable 


unacceptable 

(el. scatt. peak) 


acceptable 


FTFC 






unacceptable 

(el. scatt. peak) 


acceptable 


FTFP 


unacceptable 

(el. scatt. peak) 


poor 

(shape) 


unacceptable 

(el. scatt. peak) 


acceptable 


FTFP_BERT 


unacceptable 

(el. scatt. peak) 


acceptable 


unacceptable 

(el. scatt. peak) 


acceptable 
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